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Abstract 

Since the appearance of the classical paper of Lifshitz almost half a century ago, linear 
stability analysis of cosmological models is textbook knowledge. Until recently, however, 
little was known about the behavior of higher than linear order terms in the perturbative 
expansion. These terms become important in the weakly nonlinear regime of gravitational 
clustering, when the rms mass density contrast is only slightly smaller than unity. In 
the past, theorists showed little interest in studying this regime, and for a good reason: 
only a decade ago, it would have been an academic excercise - at scales large enough to 
probe the weakly nonlinear regime, all measures of clustering were dominated by noise. 
This is no longer the case with present data. The purpose of this talk is to provide a 
brief summary of recent advances in weakly nonlinear perturbation theory. We present 
analytical perturbative results together with results of N-body experiments, conducted to 
test their accuracy. We compare perturbative predictions with measurements from galaxy 
surveys. Such comparisons can be used to test the gravitational instability theory and to 
constrain possible deviations from Gaussian statistics in the initial mass distribution; they 
can be also used to study the nature of physical processes that govern galaxy formation 
("biasing"). We also show how future studies of velocity field statistics can provide a new 
way to determine the density parameter, fi. 
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1 Perturbation theory 



Tests of theories for the origin of the large scale structure of the universe ultimately depend on 
a comparison of model predictions with measurable quantities, derived from observations. The 
statistical measures we will discuss here are the low order N-point correlation functions. They 
have two clear advantages. First, they can be estimated from galaxy surveys with a reasonable 



degree of precision and reproducibility (see, e.g. |E| |46|, |20], |5T] , and references therein). Second 



correlation functions can be relatively easily related to dynamics (e.g. f45|). Their evolution 
can be studied by taking moments of the hydrodynamical equations of motion of an expanding 
self-gravitating pressureless fluid with zero vorticity, which is a good approximation of the 
real universe after the hydrogen recombination. Low order correlation functions can also be 
measured from N-body experiments. 

1.1 Gravitational instability 

The full description of gravitational instability is nonlinear. The density contrast 

*(x,f) = ^$)-l, (1) 

the peculiar velocity v and the gravitational potential </> are related by the Euler, Poisson and 
continuity equations. They can be combined into one expression for the density contrast, 

8 + 2H8 - \ttH 2 5 = \VtH 2 5 2 + a - 2 V8 ■ V</> + a~ 2 V a Vp[(l + 8)v a v ] . (2) 

Here x = {x a } and t are, respectively, the comoving spatial coordinates and the cosmological 
time, the dots represent time derivatives, V Q = d/dx a , a(t) is the scale factor, H = a/a is the 
Hubble parameter, and Q = 8irG(p) /3H 2 is the density parameter. The brackets (. . .) denote 
ensemble averaging and p is the mass density. Perturbation theory rests on a conjecture that 
when the deviations from homogeneity are small, the first few terms of the expansion 

5 = s 1 + 8 2 + 8 3 + ... (3) 

provide a reasonable approximation of the exact solution of eq. (0). The first, linear term is 
the well known Lifshitz |39| solution of eq. (§) with the right-hand side set to zero, 



8i = -D(£)e(x) + decaying mode , (4) 

where D(t) is the standard growing mode (see, e.g., §11 in |i5|| ) and e is a random field with 
statistical properties defined by initial conditions. The term 5 2 — ^(^ 2 ) is the solution of the 
same equation with quadratic nonlinearities included iteratively by using 8i as source terms 
(as in [45|, §18). For a vanishing cosmological constant (A = 0) and arbitrary Q, the fastest 



growing mode in this solution is [|35| 

5 2 = D 2 (t) [l(l + K)e 2 + Ve- V$ + (f - «)t^t q/3 ] , (5) 

r a(i = (|cWV 2 - V Q V^)$ ; $(x) = - / dVe(xO/47r|x - x'| . (6) 

The parameter k is a slowly varying function of time; its dependence on Q is extremely weak, 
too. In the range 0.05 < Q < 3, the ^-dependence is well approximated by 

k [m\ « n ^ 2/63 • (7) 

The higher order terms can be constructed out of the linear solution in a similar way. 
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1.2 Dynamical evolution of correlations 



The main source of difficulties in studies of the dynamics of gravitational instability is the 
nonlinearity of the equations of motion. Indeed, let us write eq. (||) for <5(x, t) and 5' = 5(x', t), 
then multiply the first equation by 8', and the second by 5. The result of adding these two 
equations and averaging is 

'( + 2Hi-mH 2 i + 2a- 2 V 2 i v = 3VlH 2 (5 2 5') + ... , (8) 

where £ = (55') and fa = (v • v') are the two-point density and velocity correlation functions, 
while the right-hand side of the equation involves 3- and 4-point correlations, which we left out 
except for one term given as an example. The full expression is not the point here; the important 
message is that we have an infinite hierarchy of equations. Mathematically, this is similar to 
the BBGKY hierarchy: in order to solve for an iV-point function, we need to know (N+l)- and 
(iV + 2)-pomt correlations. Perturbation theory may help, but not without further assumptions, 
necessary to close the infinite hierarchy of correlation functions. Here we will assume that e(x) 
is a random Gaussian field. This approach, also known as the random phase approximation, 
can be justified theoretically for a wide class of models of the early universe [@. Independently 
of whether we take "inflationary" arguments seriously or not, the random phase assumption is 
observationally testable. As we will show below, comparing model predictions with observations 
can help us decide whether the early universe was Gaussian or not. The greatest advantage of 
our assumption is simplicity. Indeed: all odd-order correlation functions of e vanish, and all 
even-order moments can be reduced to the two-point function. For example, for the first four 
moments we have ( 1 ) = ( 123 ) = 0, ( 12 ) = fa (|xi — x 2 |, t) D~ 2 , and 

(1234) = J D- 4 a(|x 1 -x 2 |,t)£ L (|x3-x 4 |,t) + cycl. (two terms) , (9) 

where (12... N) = ( £(xi)e(x 2 ) . . . e(xjy) ), while fa in eq. (§) is ( 5i (x 1; t) 5\ (x 2 , t) ), the 
linear 2-point correlation function. Using these relations and the perturbative expansion, we 
can express the 3- and 4-point functions on the right-hand side of eq. (§) in terms of fa. This 
will close the hierarchy, allowing us to calculate weakly nonlinear corrections for the 2-point 
function £. The same procedure can be applied for higher order correlation functions. This 
simple idea is at the root of all perturbative results, discussed below. 

Perturbative calculations are usually easier after applying the Fourier transform. The 
Fourier transform of £ is the power spectrum, P(k,t), where k is the comoving wavenumber. 
The linear term in the perturbative expansion for P(k, t) is 

P L {k,t) = J fa(r,t) exp(zk-r)d 3 fc . (10) 

The basic properties of linear evolution, as well as the effect of nonlinear interactions can be 
deduced from the functional form of 5\ and 5 2 . Indeed, according to eq. @), all fluctuations 
evolve without changing the shape of their density profiles, so that we can expect that the 
linear power spectrum evolves without changing its slope (the growth rate is the same for 
all scales). This is indeed the case. If Pi,(k,t) oc k n at some initial time to, than at later 
times Pi(k,t) = Pl (k,t ) D 2 (t)/ D 2 (t ) oc k n . A second characteristic property of a linear 
field is that its distribution function remains Gaussian; only its width, (5\) grows in time as 
D 2 (t). Neither of these properties remain valid when nonlinear interactions become important. 
The nonlinear terms, describing tidal interactions in eq. (|5|) redistribute power across the 
spectrum, and growth rates may be enhanced at some wavenumbers and suppressed at others: 
the logarithmic slope n(k,t) = d(log P)/d(logfc) is no longer constant, it can now change with 
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time. Moreover, the nonlinear evolution introduces deviations from Gaussian statistics. The 
lowest-order moment which can detect such deviations is the skewness, (8 3 ). The Gaussian 
distribution is symmetric about 5 = 0, and its skewness is zero. According to eq. (M), however, 
the fluctuations with positive initial density contrast 8\ grow more rapidly than those with 
5i < 0. Hence, nonlinear interactions introduce asymmetry in the distribution, and we can 
expect that with time, gravity will induce a non-zero skewness in the distribution. 

The earliest (to our knowledge) analytic results of weakly nonlinear calculations were ob- 
tained by Zeldovich et al. and Peebles [|56|, [19], ^5] . They considered power-law initial conditions, 
Pl oc k n , and showed that unless n < 4, the nonlinear terms become larger than linear at large 
scales, and the perturbation theory breaks down. A somewhat more systematic approach, go- 



ing beyond scaling relations was introduced in the eighties by Juszkiewicz et al. |33], |34| and 
Vishniac []53|| , who used Eulerian perturbation theory and the random phase approximation 
to study the first non-trivial, weakly-nonlinear corrections for second moments - P(k, t) and 
£(r, t). Peebles used the same approach to calculate the gravitationally induced skewness, 
while Fry, in a truly seminal paper, showed how to calculate a reduced correlation function of 
arbitrary order N, using the perturbative expansion, truncated at the (N — 1) -st order |23| . 

Unfortunately, N-body simulations, available in the early eighties were too crude to study 
the range of validity of perturbative methods. Moreover, the galaxy surveys at scales, large 
enough to probe the weakly nonlinear regime of clustering, were dominated by noise. This 
made further progress in analytic calculations difficult and the field became dormant for several 
years. With time, the dynamical range of the simulations improved by more than an order of 
magnitude [f)5|] . Simultaneously, the length scales in galaxy surveys, at which correlations can be 
accurately measured, increased by an order of magnitude [pOfl , making weakly nonlinear theory 
useful in their analysis. As a result of these developments, the study of the evolution of the 
mass distribution in the expanding universe has become a rich topic for analytic perturbative 
methods. The revival of enthusiasm for analytic calculations produced many new results. It is 
our purpose to briefly summarize some of these new results here. For a more detailed approach, 
see e.g. the recent lecture notes of Bouchet M and the review by Sahni & Coles 



1.3 Spatial smoothing 

The last issue we need to discuss before moving on to the next section is spatial smoothing. It 
is useful to employ a field 5, such that 5(x, t) is replaced by 5± + 8 2 + . . ., with 

8 N (x,t) = J 8 N (x.',t)w(\x-x!\/R)d 3 x' , (11) 

where w(x/R) is a spherically symmetric window function. R is the characteristic co moving 
smoothing length: w decreases rapidly for x/R ^> 1. Its Fourier representation acts as a low 
pass filter, suppressing the fuctuations with wavenumbers k 3> 1/R. The usual normalization 
condition requires / d 3 rw(r) = 1. Smoothing over a ball of radius R is also called "top hat" 
smoothing; another popular filter is a fuzzy ball with a Gaussian profile of width R. Note that 
the effects of smoothing and gravitational evolution are interchangeable only for 8i; for terms 
of second or higher order, these two effects do not commute (this is evident from eq. (|5|); see 
also |3(|). For a power-law initial spectrum, the smoothed linear variance is 

(5?) = (?l{R,t) = J w(x 1 /R)w(x 2 /R)Z L (\x 1 -x 2 \)d 3 x 1 d 3 x 2 oc R' {n+3) , (12) 

so the condition n > —3 ensures that the clustering process is hierarchical, and proceeds from 
small to large scales, (this constraint, together with the Zeldovich-Peebles condition n < 4, 
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broadly determines the range of "reasonable" values of n; a mix of observational arguments and 



theoretical prejudce shrinks this range to —2 <n <1, see e.g. [46]). The smoothing procedure 
is needed for at least two reasons. First, at any stage of the evolution, in order to be able to use 
the perturbation theory, we need to smooth-out the strong nonlinearities. In all hierarchical 
models, this can be done by any a low pass filter of width R > R nh where R n i(t) is the current 
nonlinear scale, defined by the condition o~ 2 L (R n i,t) = 1. The second reason for introducing 
the smoothing is practical: all observable quantities, derived from galaxy surveys or numerical 
experiments, involve spatial smoothing. Hereafter, we will deal with smoothed fields almost 
without exceptions (which will be explicitly stated), so to simplify notation, we will omit the 
bars over 5 and 5n- Generally, we will also assume Gaussian initial conditions (again, all 
exceptions will be properly described). 



2 New results for the second moment; the previrialization conjec- 
ture 

The perturbative series (§) can be used to expand the second moment of the density field in 
powers of the linear variance, o\. Squaring 5 and averaging, we get 

a 2 = (5 2 ) = (S 2 ) + (6 2 2 ) + 2(S 1 S 3 ) + . . . = a 2 L + (/ 22 + 7 13 ) a* + 0(a 6 L ). (13) 

The second and third term in this expansion are the lowest order nonlinear corrections to linear 
theory. These terms are of order aj^ because O(o~\) terms vanish under the random phase 
assumption. The quantities I22 and /13 can be determined from perturbation theory. They 



were recently calculated for power-law initial conditions by Scoccimarro & Frieman |50|] and 



by Lokas et al. [41]. Scoccimarro & Frieman explored the evolution of an unsmoothed density 
field. Lokas et al. studied spatially smoothed fields, and found that nonlinear tidal interactions 
redistribute the power in the spectrum in such a way, that for n > — 1, the correction 7 2 2 + 113 
is negative. The growth rate of density fluctuations is inhibited: a is smaller than the linear 
prediction <jl. The magnitude of this suppression increases with n, or with the relative amount 
of small scale power in the initial spectrum. For n < — 1, there is an opposite effect - the 
fluctuation growth rate is enhanced (a > <Jl). The transition occurs at n = —1, when the 
correction is close to zero and a ~ a L . Lokas et al. compared these perturbative results with 
a series of N-body simulations and found a good agreement in a wide dynamical range, up to 
a ~ 1. These results are in qualitative agreement with the so-called previrialization conjecture 



of Davis & Peebles ||18|| . According to Davis & Peebles, tidally induced nonradial motions 
of small scale inhomogeneities within collapsing mass concentrations should resist gravity and 
slow down the collapse. Until recently the size of the previrialization effect was still under 
discussion, with some N-body studies reproducing the effect and others claiming its absence. 
For example, Peebles [57] found strong support for the previrialization effect, while Evrard Sz 



Crone |21| reached opposite conclusions. This apparent controversy can be resolved, using the 
perturbative results we just discussed. Indeed, Peebles in his simulations assumed n = 0, while 
Evrard & Crone assumed n = —1: a spectral index, for which nonlinear effects are minimal. 



Jain & Bertschinger []32[ used perturbation theory and N-body experiments to investigate 
the nonlinear transfer of power in a CDM spectrum. Their results are consistent with Lokas et 
al. and Scoccimarro & Frieman. Finally, we must mention the important work of Makino et al. 



4*2*|, who succeeded in reducing the dimensionality of the mode coupling integrals of Vishniac 



and Juszkiewicz et al. J53], 34|. This significantly simplified the derivation of all of the new 



results described in this section. 
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3 The third moment and deviations from Gaussian behavior 



As we saw in Section 1, nonlinear gravitational evolution drives the field 5 away from the initial 
Gaussian distribution. Deviations from Gaussianity can be measured by cumulants, also called 
"reduced" or "connected" moments. A cumulant of order N is 



M N = d N \n(e xd )/dx N , evaluated at x = . (14) 

For the first four cumulants, this expression generates the mean, M± = (5) = 0, the variance, 
M 2 = (S 2 ) = a 2 , the skewness, M 3 = (<5 3 ), and the kurtosis, M 4 = (5 A ) — 3<r 4 . For a zero-mean 
Gaussian distribution, all cumulants vanish except M 2 . The skewness measures the asymmetry 
of the distribution, while the kurtosis measures the flattening of the tails relative to a Gaussian. 
Using perturbation theory, Fry [^] showed that in an Einstein-de Sitter universe, Mn scale 



like a 2N 2 . Hence, to lowest non- vanishing order in perturbation theory, the ratio 

S N = M N /a 2N - 2 (15) 



is time- independent. Juszkiewicz et al. [[36j have recently used this property and the Edge- 
worth series to describe the weakly nonlinear evolution of the probability distribution function 
(PDF) of the density field. The Edgeworth expansion is a standard statistical tool, used to 
represent quasi-gaussian probability distribution functions (PDFs). Juszkiewicz et al. obtained 
an expansion in powers a, with the first few terms, given by 

p(u) = {l + ± S 3 H 3 (v)a + {± S 4 H 4 (v) + §S 2 H 6 (v)}a 2 } g(u) + O (g(v)a*) , (16) 

where v = 5/cr, g = (27r) -1 / 2 exp(— z/ 2 /2), and Hj is a Hermite polynomial of degree J. A 
similar expansion was obtained independently by Bernardeau & Kofman 0. The Edgeworth 
expansion is particuraly useful when only a few low-order reduced moments are known; it can 
be also used to estimate low-order Sn -s from the shape of empirically measured PDF near the 
peak (in realistic situations, the tails, and/or high-iV moments are poorly known anyway). An 
alternative, Laplace-transform representation for the PDF, which does require the knowledge 
of all Sn -s to iV = oo, was proposed by Bernardeau 0. 

3.1 Skewness induced by gravity 

The acual values of the parameters Sjy can be also calculated perturbatively. The results 
generally depend on the shape of the initial power spectrum and the shape of the window 
function. Let us begin with the skewness, responsible for the lowest order nonlinear correction 



to the Gaussian PDF in eq. (|T6|) . Upon raising both sides of equation (|3]) to the third power 
and averaging, we get 

(5 3 ) = S 3 a 4 = (6*) + (36152) + ... (17) 

In models, which satisfy the random phase assumption, the lowest order non-vanishing term 
in this expansion is O (cr) 4 . The scaling with a is completely different in models with non- 
gaussian initial conditions and an intrinsic skewness (5f) ^ 0: the lowest order term is oc a 3 . 



The skewness parameter, S3 in eq. ( \F7\) was first calculated by Peebles [45], assuming f2 = 1 and 
no smoothing. The latter constraint prevents any comparison of this result with observations 
or N-body simulations. To deal with this problem, we calculated S3 for a smoothed density 
field with a power-law initial spectrum and for two different kinds of spatial filters: a gaussian 



and a top hat p5| . The result is particularly simple in the latter case, 



S 3 = f - (n + 3) + 0(a 2 ), (18) 
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where the allowed range of power indexes is —3 < n < 1 (outside of this range S3 diverges). 
Note that here, as well as in Section 2, we quote only results for the Einstein-de Sitter case. 
All these results, however, are also true for a much wider class of models with arbitrary Q and 
A (this is an approximate statement, but the approximation is good enough for all practical 
purposes). For example, for < A/3H 2 < 1 — Q and 0.05 < Q < 3 , the term '34/7' in eq. 
( |T8[ ) should be replaced by 

4 + 4k(0) = f + f (H- 03 - 1) . (19) 

The corrections, introduced by deviations from the Einstein-de Sitter model were first derived 
by Bouchet et al. || for arbitrary Q, and by Bernardeau || and Bouchet et al. for arbitrary 
A ||10|| . These results were recently confirmed by Catelan et al. |14|], who used a significantly 



different method of derivation. For our purposes here, the most important implication of all 
this is that Gaussian initial conditions and weakly nonlinear gravitational instability generate 
the scaling relation 

M N oc a 2N - 2 (?<1), (20) 

and Sn, the ratio of the above moments, up to terms of order a 2 in perturbation theory is 
time-, Q- and A-independent; it is determined only by the shape of the initial power spectrum. 
Bernardeau generalized eq. ([18]) for a smoothing scale-dependent slope, n = n(R). He 
showed, that this expression for Ss(R) remains valid when n is replaced by the effective local 
slope, 

n(R) = -d(loga£0R))/d(logfl)-3. (21) 

As we have seen, to lowest non- vanishing order in perturbation theory, S3 is time independent. 
We can expect that S 3 (i?) will be scale-independent as well in the range of smoothing scales, 
for which n(R) is constant. A completely different picture arises when the primeval PDF is 
significantly skew. The lowest order non-vanishing contribution to skewness gives 

S 3 (R) oc l/a L oc R^ n+z)/2 D{t)~ l , (22) 

so in all hierarchical (n > —3) models we can expect a strong scale- and time-dependence, 
with S3 "blowing up" at large scales and at high redshift. As pointed out by several authors, 
this property can be used to distinguish observationally Gaussian initial conditions from their 
non-gaussian alternatives [T7|, |35], [EJ. Chodorowski & Bouchet [[UJ find that the distinction 
can be made even sharper by studying the scaling properties of S4 (see also Fry & Scherrer 
251V 



Using methods similar to those discussed above, Catelan & Moscardini [15], Lokas et al. 
i0" |, and Bernardeau [|2[ calculated the parameter S4 for power-law spectra and top hat as well 
as Gaussian windows. The last two papers also provide closed-form expressions for S4 = S^n). 
Bernardeau || solved the problem of calculating Sn -s of arbitrarily high order for top hat 
smoothing. 

Most of the above perturbative results for S 3 and S 4 were tested against N-body simulations 
[[I], 0, [35], |38|, |40], showing remarkable agreement up to a w 1. On a comoving scale 



.R, the Edgeworth expansion agrees with the N-body simulations within ~ l/a(R) standard 
deviations from zero J3BJ. The Laplace transform perturbative PDF agrees with the N-body 
PDF even when o is slightly greater than unity j5]. 

3.2 Velocity statistics 

Unlike the reduced moments of the density field, the cumulants of the expansion scalar, 9 = 
V ■ \/H are very sensitive to the value of fl For example, the skewness parameter is given by 
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the expression |J 

(9 3 )/(9 2 ) 2 = -ft- - 6 [M-( n + 3)] . (23) 

In principle, this can provide a new method of determining Q from the velocity data alone 
(rather than from the comparison with the density field, as in more conventional methods). 
The existing velocity surveys, however, are by far too shallow to allow accurate measurements 
of the variance and skewness of 9. We have to wait for better data. Weakly nonlinear theory 
can also be used to study the dynamics of relative motions in pairs of galaxies, which induce 



the anisotropy of the redshift space galaxy correlations [57]. This can improve the accuracy of 



fl determinations from presently available redshift surveys. 
3.3 Lagrangian methods 

So far, we have focused on Eulerian perturbation theory. An interesting alternative is the 
Lagrangian approach. It follows particle trajectories rather than the evolution of 8(pc, t). The 
local density is <5 = 1/ J — 1, where J is the Jacobian of the transformation from initial particle 



positions to positions at time t. The famous Zeldovich approximation |57J is a linear order 
Lagrange solution. In its original form, it can not be used for statistical studies because 
it predicts infinite cumulants for 5. As pointed out by Grinstein & Wise PD| , the problem 
can be at least partially fixed by expanding 1/ J — 1 in powers of spatial derivatives of the 
linear order Eulerian velocity field. The cumulants, generated by such an expansion are finite. 
Unfortunately, however, they disagree with the proper Eulerian perturbation theory and with 
all N-body experiments |~5], |35], |36| . The Zeldovich approximation fails particularly badly when 



applied to the velocity divergence, 9 S !5J]. Two lessons follow from this example. The first 
is that dynamical estimates of Q, derived from reconstruction schemes for the initial PDF of 
5 or 9, which rely entirely on the Zeldovich approximation (see e.g. ||44||), should be treated 
with caution. The second lesson is positive: the Zeldovich approach is more powerful than it 
seems. For example, under this approximation, k = in eq. ([19]) and therefore for a top hat 
S3 = 4 — {n + 3). Depending on n, this causes about 20% to 50% underestimate of the true S3, 
derived from second order Eulerian theory. Note however, that we used first order Lagrangian 
theory to derive the Zeldovich prediction for 5*3. To make a fair comparison, we should confront 
this result with its first order Eulerian counterpart: the linear prediction S3 = 0, which is 
much further away from the true value! What we have just seen, is an illustration of a more 
general phenomenon: at any fixed order in perturbation theory, the Lagrangian approach is 
likely to remain valid for stronger density contrasts than the Eulerian approach. This happens 
because the requirement of the smallness of the (Lagrangian) particle displacements is less 
restrictive than the requirement of the smallness of \S\. This idea, which must have motivated 
Zeldovich, remains valid at higher orders in perturbation theory. The idea to extend the 
Lagrangian perturbation theory beyond the linear order was introduced by Moutarde et al. 



[]43| and generalized by Bouchet et al. || and Buchert & Ehlers ][H|. Lagrangian coordinates 
have been used to derive second- and third-order perturbative solutions in the Einstein-de Sitter 
case [12], O] , as well as in more general cosmological models 0, [10[ , and to study distortions of 



correlations in redshift space |L0|, |31fl - We are convinced that the full potential of this approach 



to perturbation theory is far from being exhausted. 
3.4 Convergence studies 

As we have seen, perturbative expansions for (5 2 ) and (<5 3 ), truncated at the lowest-order or 
next to the lowest-order non-vanishing contribution, are in remarkable agreement with N-body 
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simulations in a wide dynamical range, up to a ~ 1. The precise range of validity of perturba- 
tion theory is still under investigation. Apart from comparisons with N-body experiments, the 
convergence can be also studied within the perturbation theory itself, by comparing contribu- 
tions coming from terms of increasing order. Such studies may answer the question: why does 
perturbation theory work so well? 



Scoccimarro & Frieman |5(J obtained closed-form expressions for I 2 = I22 + I13 m the case 
of no smoothing (R = 0). Although individually, I22 and J13 suffer from "infrared" (k — > 0) 
divergencies, their sum is finite and approximately n-independent; for —2 < n < 2 , the 



correction term is 41, 50 



I 2 « 1.82. (24) 

The divergent terms in I22 and /13 cancel each other. The probable source of this phenomenon 
is the Galilean invariance of the equations of motion | |50| | . 

In order to make a quantitative comparison with N-body results possible, Lokas et al. 



included the effects of smoothing in their calculations. Introduction of a finite smoothing length 
[R 7^ 0) generates "ultraviolet" {k — > 00) divergences at n > — 1; 1% diverges unless -Pl(&) has 



a shortwave cutoff at k = k c . For k c R ^ 1, the correction term is |41[ 

h « -{0.6k c R) n+1 . (25) 

We can expect that perturbation theory breaks down when <jl reaches the value, at which the 
first two terms in the expansion ( |13D become comparable, |/2|o"i = c£- Therefore, equation 
fl25|) suggests that for a given k c R, the maximal value of gl, determining the range of validity 
of perturbation theory, should gradually decrease with increasing n. Comparisons with N-body 
experiments show that this is indeed the case Note that a shortwave cutoff naturally 



appears in N-body experiments; k~ x corresponds to the mean interparticle separation, or the 
Nyquist scale. 

Scoccimarro & Frieman [^0] have also calculated the 0(a 2 ) term in the expansion for S 3 



for the unsmoothed density field. They found that 0(a1) terms in the expansions for S3 do 
not "overtake" the lowest-order non- vanishing terms until <jl ~ 0.5, in qualitative agreement 
with the results of direct N-body tests of perturbation theory, discussed above. A quantitative 
comparison, however, will be possible only after smoothing effects are included (Scoccimarro & 
Frieman, in preparation). 



4 Observations 

Before comparing theoretical predictions to observations, let us discuss two possible sources of 
confusion. 

1. Biasing. If galaxies form with different efficiency in different environments, then the 
galaxy distribution may paint a biased picture of the underlying mass distribution. Let us 
suppose that the smoothed galaxy field 5 g is a local, but not necessarily linear function of the 
smoothed mass density field: 5 s (x) = /[5(x)], where / is an arbitrary function. Then, using a 
Taylor expansion for \5\ 1, it is possible to prove, as Fry & Gaztanaga did, that local 
bias preserves the form of the scaling relations (^), while changing the constants S^. For 



example, 5*3 changes to |24], |36| 



S 3g = (5 s g ) I (5 2 2 ) = S 3 /h + 3b 2 /bi , (26) 

where bj = d J f/d5 J , evaluated at 5 = (note that the linear bias model with all bj = for 
J > 1 must break down when b > 1 and 5 < —1/6, since it predicts negative S g ). 
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2. Redshift distortions. Redshift surveys use radial velocities of galaxies instead of 
their true radial coordinates. As a result, galaxy peculiar motions distort the true spatial 
distribution and the two-point correlation function, £ (cf. §20 in and references therein). 



Bouchet et al. |L0| and Hivon et al. have used weakly nonlinear perturbation theory and N- 
body simulations to study the effects of redshift space distortions on the three-point correlation 
function, (, and on 5*3. They showed that although (, like £, is affected by the redshift space 
distortion, all appreciable, fi-dependent effects cancel out for the moment ratio S3. In the 
weakly- nonlinear domain, when o < 1, and for spectral slopes —3 < n < — 1, the differences 
between redshift-space 5*3 and its true value are smaller than 10%, whatever the value of Q. 
Hence, the comparison between theoretical predictions and observations is not obscured by 
redshift space distortions provided we use the moment ratio S3 rather than reduced moments 
themselves. 



The parameters 63 and S4 have been estimated from the CfA and SSRS surveys |pq| , the 
IRAS 1.2 Jy sample |§, as well as for vaious angular catalogues |27], |52|. All these measurements 
show beautiful agreement with the predicted scaling relation (^). The dependence (|26| ) of S 3g 
on the shape of the bias function may explain why the value measured by Bouchet et al. |J 
from the top-hat smoothed, IRAS redshift survey is rather low, S^ AS ~ 1.5. The galaxies 
observable in the infrared are underrepresented in rich clusters relative to optical galaxies, 
so the "bias function" that applies to IRAS galaxies may have negative curvature (62 < 0), 
pushing S3 9 below the value of S3 for the mass. If cluster avoidance is indeed the cause of 
the low value of S;^ 45 , then one would expect that optically selected surveys would yield 
higher values. The values of Ss g , obtained recently by Gaztanaga |27j from the APM survey 
suggest that this is exactly the case. The value of 633, deprojected from the two-dimensional 
APM data onto three-dimensional space is 3.16 ± 0.14 , assuming a top-hat filter with a radius 
R > 7/i~ 1 Mpc (here h is the Hubble constant in units 100kms~ 1 Mpc _1 ; the latter condition 
ensures cr g < 1). Substituting n = —1, appropriate for the slope, seen in the APM survey 
on these scales, into equation (|T8|) yields S3 = 2.9 . Hence, for sufficiently large scales, the 



APM results are in perfect agreement with second order perturbation theory: S^ PM = S3 . 
This agreement provides support for three theoretical assumptions: first, that gravitational 
instability is responsible for the growth of clustering; second, that the initial conditions were 
Gausssian; and third, that the APM galaxies trace the mass, i.e., in eq. ( p6|) b APM « 1 and 
b2 PM ~ (of course, other combinations of b\ and bi also fit the APM data, but they require 
a coincidental compensation of the linear, \jb\ factor by the 62 term arising from nonlinear 



bias). This is bad news non-local biasing models like the model of Bower et al. JTTJ] , invoked 
in a desperate attempt to explain the mismatch of the observed power spectrum with standard 
CDM. Indeed, as demonstrated by Frieman & Gaztanaga H22fl , such models break the scaling 
relation (BDT) and are incompatible with observations. Most recently, Gaztanaga & Mahonen 



28| have used the APM data to severely constrain the so-called global texture model. The 
source of problems for this model is its prediction that Sj parameters must 'blow up' at large 
scales, like S3 in our example in eq. fl22"|) . There is no trace of such behaviour in the APM 
galaxy distribution. Since the already available data can be used to constrain non-local biasing 
schemes and deviations from Gaussian distribution, we can expect that in the new generations 
of catalogues (like the SDSS and the 2dF, described in these proceedings) will finally tell us 
whether all we are looking at is just a Gaussian random process modified by plain gravity. 



Acknowledgements. We are grateful to the Organizers, Chantal Balkowski, Sophie Mau- 
rogordato, and J. Tran Thanh Van, for making this scientifically useful meeting happen; to 
Mother Nature for providing good skiing conditions, and last but not least, to Chantal and 



9 



Sophie as Editors, for their incredible patience with the authors, who broke several submission 
deadlines for this manuscript. This research was supported by the Polish Government KBN 
grants No. 2-P304-01607, 2-P304-01707, and by an NSF grant No. PHY94-07194. RJ thanks 
Alain Omont for his hospitality at Institut d'Astrophysique de Paris, CNRS, where part of the 
work reviewed here was done. 

References 

[1] Baugh, CM., Gaztahaga, E., Efstathiou, G., 1994, Mon. Not. R. astr. Soc. 274, 1049. 
[2] Bernardeau, F., 1994, Ap. J. 433, 1. 

[3] Bernardeau, F., Juszkiewicz, R., Dekel, A., Bouchet, F.R., 1995, Mon. Not. R. astr. Soc. 
274, 20. 

[4] Bernardeau, F., & Kofman, L., 1995, Ap. J. 443, 479. 
[5] Bernardeau, F., 1994, Astr. Astrophys. 291, 697. 

[6] Bouchet, F. R. 1995, in: Proc. of the "Enrico Fermi" School "Dark Matter in the Universe" 
held in Varenna (Societa Italiana di Fisica). 

[7] Bouchet, F.R., Hernquist, L., 1992, Ap. J. 400, 25. 

[8] Bouchet, F.R., Juszkiewicz, R., Colombi, S., Pellat, R., 1992, Ap. J. (Letters) 394, L5. 

[9] Bouchet, F.R., Strauss, M., Davis, M., Fisher, K.B., Yahil, A., & Huchra, J. P., 1993, Ap. 
J. 417, 36. 

[10] Bouchet, F.R., Colombi, S., Hivon, E., & Juszkiewicz, R., 1995, Astr. Astrophys. 296, 
575. 

[11] Bower, R.G., Coles, P., Frenk, C.S., & White, S.D.M., 1993, Ap. J. 405, 403. 

[12] Buchert, T., 1994, Mon. Not. R. astr. Soc. 267, 811. 

[13] Buchert, T., & Ehlers, J., 1993, Mon. Not. R. astr. Soc. 264, 375. 

[14] Catelan, P., Lucchin, F., Mattarrese, S., Moscardini, L., 1995, Mon. Not. R. astr. Soc. 
276, 39. 

[15] Catelan, P., & Moscardini, L., 1994, Ap. J. 426, 14. 

[16] Chodorowski, M., & Bouchet, F.R., 1996, Mon. Not. R. astr. Soc, in press. 

[17] Coles, P., Frenk, C.S., 1991, Mon. Not. R. astr. Soc. 253, 727. 

[18] Davis, M., & Peebles, P.J.E., 1977, Astrophys. J. Suppl. Ser. 34, 425. 

[19] Doroshkevich, A.G., Zeldovich, Ya. B., 1975 Astrophys. Space. Sci. 35, 55. 

[20] Efstathiou, G., 1996, in: Les Houches Lectures, Session LX, ed. R. Schaeffer, Elsevier 
Science Publishers (Netherlands). 

[21] Evrard, A.E., & Crone, M.M., 1992, Ap. J. (Letters) 394, LI. 

[22] Frieman, J., & Gaztahaga, E., Ap. J. 425, 392. 

[23] Fry, J.N., 1984, Ap. J. 279, 499. 

10 



[24] Fry, J.N., & Gaztanaga, E., 1993, Ap. J. 413, 447. 

[25] Fry, J.N., & Scherrer, R.J. 1994, Ap. J. 279, 36. 

[26] Gaztanaga, E., 1992, Ap. J. (Letters) 398, L17. 

[27] Gaztanaga, E., 1994, Mon. Not. R. astr. Soc. 268, 913. 

[28] Gaztanaga, E., & Mahonen, P., 1996, Ap. J. (Letters), in press 

[29] Grinstein, B., & Wise, M.B., 1987, Ap. J. 314, 448. 

[30] Goroff, M.H., Grinstein, B., Rey, S.-J., & Wise, M.B., 1986, Ap. J. 311, 6. 

[31] Hivon, E., Bouchet, F.R., Colombi, S., & Juszkiewicz, R., 1995, Astr. Astrophys. 298, 
643. 

[32] Jain, B., Bertschinger, E., 1994, Ap. J. 431, 495. 

[33] Juszkiewicz, R., 1981, Mon. Not. R. astr. Soc. 197, 931. 

[34] Juszkiewicz, R., Sonoda, D.H., & Barrow, J.D., 1984, Mon. Not. R. astr. Soc. 209, 139. 

[35] Juszkiewicz, R., Bouchet, F.R., & Colombi, S., 1993, Ap. J. (Letters) 412, L9. 

[36] Juszkiewicz, R., Weinberg, D.H., Amsterdamski, P., Chodorowski, M., & Bouchet, F., 
1995, Ap. J. 442, 39. 

[37] Juszkiewicz, R., Fisher, K.B., & Szapudi, I., 1996, Ap. J., in press. 

[38] Lahav, O., Itoh, M., Inagaki, S., & Suto, Y., 1993, Ap. J. 402, 387. 

[39] Lifshitz, E.M., 1946, J. Phys. U.S.S.R. 10, 116, (Zh.E.T.F. 16, 587). 

[40] Lokas, E.-L., Juszkiewicz, R., Weinberg, D.H., & Bouchet, F.R., 1995, Mon. Not. R. astr. 
Soc. 274, 730. 

[41] Lokas, E.-L., Juszkiewicz, R., Bouchet, F.R., & Hivon, E., 1996, Ap. J., in press. 
[42] Makino N., Sasaki M., Suto, Y., 1992, Phys. Rev. D46, 585. 
[43] Moutarde, F., Alimi, J.-M., Bouchet, F.R., & Pellat, R., 1991, Ap. J. 382, 377. 
[44] Nusser, A., & Dekel, A., 1993, Ap. J. 405, 437. 

[45] Peebles, P.J.E., 1980, The Large-Scale Structure of the Universe (Princeton University 
Press: Princeton). 

[46] Peebles, P.J.E., 1993, Principles of Physical Cosmology (Princeton University Press: 
Princeton) . 

[47] Peebles, P.J.E., 1990, Ap. J. 365, 27. 

[48] Sahni, V., & Coles, P., 1995, Physics Reports, in press (preprint SISSA |astro-ph/9505005l) . 

[49] Silk, J., & Juszkiewicz, R., 1991, Nature 353, 386. 

[50] Scoccimarro, R., &; Frieman, J., 1995, Fermilab preprint. 

[51] Strauss, M.A., & Willick, J.A., 1995, Physics reports 261, 271. 

[52] Szapudi, I., Dalton, G., Efstathiou, G., & Szalay, A.S., 1995, preprint. 

11 



[53] Vishniac, E., 1983, Mon. Not. R. astr. Soc. 203, 345. 

[54] Weinberg, D.H., & Cole, S., 1992, Mon. Not. R. astr. Soc. 259, 652. 

[55] White, S.D.M., 1996, in: Les Houches Lectures, Session LX, ed. R. Schaeffer (Elsevier 
Science Publishers: Netherlands). 

[56] Zeldovich, Ya.B., 1965, Adv. Asron. Astrophys. 3, 241. 
[57] Zeldovich, Ya.B., 1970, Astr. Astrophys. 5, 84. 



12 



